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We calculate the coefficients in the chiral Lagrangian approximately from QCD based on a 
previous study of deriving the chiral Lagrangian from the first principles of QCD in which the 
chiral Lagrangian coefficients are defined in terms of certain Green's functions in QCD. We first 
show that, in the large- N c limit, the anomaly part contributions to the coefficients are exactly 
cancelled by certain terms in the normal part contributions, and the final results of the coefficients 
only concern the remaining normal part contributions depending on QCD interactions. We then 
do the calculation in a simple approach with the approximations of taking the large- iV c limit, the 
leading order in dynamical perturbation theory, and the improved ladder approximation, thereby 
the relevant Green's functions are expressed in terms of the quark self energy S(p 2 ). By solving the 
Schwinger-Dyson equation for S(p 2 ), we obtain the approximate QCD predicted coefficients and 
quark condensate which are consistent with the experimental values. 

PACS number(s): 12.39.Fe, 11.30.Rd, 12.38.Aw, 12.38.Lg 



I. INTRODUCTION 



Because of its nonperturbative nature, studying low energy hadron physics in QCD is a long standing difficult 
problem. For low lying pseudoscalar mesons, a widely used approach is the theory of the effective chiral Lagrangian 
based on the consideration of the global symmetry of the system and the momentum expansion without dealing with 
the nonperturbative dynamics of QCD Q In the chiral Lagrangian approach, the coefficients in the Lagrangian 
are all unknown phcnomcnological parameters which should be determined by experimental inputs. The number 
of unknown parameters increases rapidly with the increase of the precision in the momentum expansion. Therefore 
studying the relation between the chiral Langrangian and the fundamental principles of QCD will not only be the- 
oretically interesting for a deeper understanding of the chiral Lagrangian, but will also be helpful for reducing the 
number of unknown parameters and increasing the predictive power of the chiral Lagrangian. 

In a previous paper, Ref. certain techniques were developed, with which the chiral Lagrangian was formally 
derived from the first principles of QCD without taking approximations. The chiral Lagrangian coefficients are 
contributed both by the anomaly part (from the quark functional measure) and the normal part (from the QCD 
Lagrangian). In Ref. Q, all the chiral Lagrangian coefficients contributed from the normal part of the theory are 
expressed in terms of certain Green's functions in QCD, which can be regarded as exact QCD definitions of the chiral 
Lagrangian coefficients. After expanding the effective action in powers of the rotated sources (momentum expansion) , 
the effective action , up to 0(p 4 ), contributed from the normal part is of the form Jl 
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where fl is related to the nonlinearly realized meson field U by U — f2 2 ; sq, pn, vn, and an are, respectively, the 
external scalar, pseudoscalar, vector, and axial-vector sources rotated by f2; and the ZCs are terms with different Lorentz 
structures in the relevant QCD Green's functions. The obtained expressions for the chiral Lagrangian coefficients are 
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Together with the anomaly part contributions, the complete coefficients are given by 
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where the superscripts (anom) and (norm) denote the anomaly part and normal part contributions, respectively. 

In the literature, the anomaly part contributions are usually calculated by means of the heat kernel regularization 
technique However, this technique is difficult to implement in the calculation of the normal part contributions 
which contain complicated functions of the momentum, say the quark self-energy £(p 2 ), reflecting nonperturbative 
QCD dynamics (which are even unspecified in the analytical part of the calculation). In order to treat the anomaly part 
and the normal part contributions on equal footing, a certain new regularization technique feasible for the calculations 
of both parts should be developed. In this paper, we use the generalized Schwinger proper time regularization technique 
developed in Ref. || to regularize the system, which keeps the local chiral symmetry at every step in the calculation, 
and can be applied to the calculations of both the contributions from the anomaly part and from the normal part. 
Thus, in this paper, the contributions from the anomaly part and the normal part are calculated by means of the 
same technique. As the first conclusion of this unified treatment, we show that the anomaly contributions to the 
chiral Lagrangian coefficients given in Ref. W, which are independent of QCD interactions, will actually be cancelled 
by certain terms in the normal part contributions, and the final expressions for the coefficients concern only the 
remaining terms from the normal part contributions related to QCD interactions. It should be so since the coefficients 
indicate meson interactions which should be residual interactions between quarks and gluons, and thus should depend 
on QCD interactions. These contributions have not been carefully calculated in the literature. It has been shown 
in Ref. |3| that in the approximations of large-iV c limit, leading order in dynamical perturbation, and improved 
ladder approximation, the formula for Fq in Eqs.(||) reduces to the well-known Pagels-Stokar formula Q in which all 
dynamical effects from QCD are represented by the quark self-energy S(p 2 ) in the formula. In this paper, we take 
the same approximations to calculate the chiral Lagrangian coefficients (the relevant QCD Green's functions) as an 
illustration of the main feature of how QCD predicts the chiral Lagrangian coefficients. Similar to the case of the 
Pagels-Stokar formula, the relevant QCD Green's functions can all be expressed as functions of the quark self-energy 
E(p 2 ). By solving the Schwinger-Dyson equation, we obtain £(p 2 ), and thus the approximate QCD predicted values 
of the coefficients. We shall see that the obtained coefficients L\, ■ ■ ■ , L w and quark condensate are consistent with 
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the experimental values. The calculation is checked by the absence of divergences in the large- N c limit as it should 
be since the divergent meson- loop contributions are of next-to-the-leading order in the 1/N C expansion. Although 
the present approximation is rather crude, it reveals the main feature of QCD predictions for the chiral Lagrangian 
coefficients. 

This paper is organized as follows: In Sec. II, we calculate the anomaly part contributions to the Ofp 4 ) coefficients 
using the Schwinger proper time regularization technique, and the results coincide with those in Ref. |4| in the chiral 
limit. Then, in Sec. Ill, we apply the same technique to the normal part, and show generally that, in the large- N c 
limit, the anomaly part contributions to the chiral Lagrangian coefficients are exactly cancelled by the contributions 
from a piece in the normal part independent of the quark self-energy, and the contributions from the remaining 
piece in the normal part depending on the quark self-energy play the real role in the chiral Lagrangian coefficients. 
Specific approximations in the calculation of the normal part contributions and the formulae for the complete chiral 
Lagrangian coefficients in terms of the quark self-energy are given in Sec. IV. In Sec. V, we present the numerical 
calculations of the quark self-energy and the obtained values of the chiral Lagrangian coefficients. Section VI is a 
concluding remark. 



II. ON THE CONTRIBUTIONS FROM THE ANOMALY PART 



In order to see the relation between the anomaly part and the normal part contributions to the chiral Lagrangian 
coefficients, we present here the calculation of the anomaly part contributions by means of the Schwinger proper 
time regularization. We shall see that the obtained results exactly coincide with those obtained from the heat kernel 
technique (J] . Our present approach is different from that in Ref. Q in the sense that the constant constituent quark 
mass Mq is not put in by hand as is done in Ref. Q but is naturally included in the normal part solution through the 
dynamical quark mass reflecting chiral symmetry breaking. Therefore our result of the anomaly part contribution is 
to compare with that in Ref. [Q in the chiral limit. 

In the Schwinger proper time regularization, the anomaly part does not contribute to the coefficients of the 0(p 2 ) 
terms in the case corresponding to the result with Mq = in Ref. HQ. Therefore we are only going to calculate the 
anomaly contribution to the coefficients of the 0(p 4 ) terms. 

The anomaly term in the path-integral formalism is 

^anom) _ _ . x anomaly terms = _ jjV c [Tr ln(ip + J) - Tr ln(ij0 + J n )] 

= iN c [Tr\n(ip + Jq) + ft — independent term]. (5) 

The O-independent term is independent of the U field, so that it is irrelevant to the chiral Lagrangian coefficients. 
We shall only evaluate the ^-dependent term in Eq.(^|). To have a unified parametrization, we can parametrize the 
anomaly contributed effective action similar to that in Eq.(|l]), i.e., 
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The ^-dependent term in Eq.(|^) suffers from ultraviolet divergence, and we take the Schwinger proper time regu- 
larization with an ultraviolet cutoff parameter A to regularize it. To apply this regularization, we first work in the 
Euclidean space-time, and analytically continue the results to the Minkowskian space-time after the evaluation. The 
main procedure of evaluating the general functional determinant including the quark self-energy S is described in 
APPENDIX A. In the case of S^ wm \ there is no ^-dependence in Eq.(||). However, for regularizing the infrared 



x If one takes a momentum cutoff A to regularize the divergent integrals as was done in Ref. Q before putting in the constituent 
quark mass Mq, the 0(p 2 ) coefficient Fq will be proportional to A 2 (cf. Ref. Q). As has been pointed out in Ref. || this term 
is exactly cancelled by a corresponding term in the normal part contribution [cf. Eq.(74) in Ref. ||]. 
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divergence, we should replace the £ in Eqs.(Al) and (A3) by an infrared cutoff parameter k. The momentum integra- 
tion in Eq.(A3) can be explicitly carried out with a lengthy but elementary calculation. After expanding in powers of 

the external sources, we can identify the expressions for /C^ anom \ • • ■ , K.^ 1 by comparing with the form of Eq.(||), 
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Comparing with the standard form of momentum expansion to identify the 0(p ) chiral Lagrangian coefficients, 
we obtain the anomaly contribution to these coefficients 
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These exactly coincide with the results with Mq = in Ref. [Q. Note that the final expressions of the coefficients 
Li, • • • , L10 are independent of the infrared cutoff parameter ft and the ultraviolet cutoff A although these cutoff 
parameters appear in IC^ nom \ • • • , , while Hi and Hi depend on the cutoff parameters. This implies that Hi 

and Hi are not measurable quantities. With N c = 3, the values of the coefficients are (in units of 1CP 3 ) 



Li = 0.79, 
L 7 = 0.26, 



L 2 = 1.58, 
L 8 = -0.79, 



L 3 = -3.17, 
Lg — 6.33, 



L4 — L 5 — L 6 — 0, 
L w = -3.17. 
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These are to be compared with the experimental values (in units of 10 



Li = 0.9 ±0.3, 
L 7 = -0.4 ±0.15, 



L 2 = 1.7 ±0.7, 
L H = 1.1 ±0.3, 



L 3 = -4.4 ±2.5, 
L 9 = 7.4 ±0.7, 



L 4 = 0±0.5, L 5 = 2.2±0.5, 
L w = -6.0 ±0.7. 



0±0.3, 
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The numbers in Eqs.(|^) are close to the experimental results of Eqs.(|l0|) except Is? and Lg are of wrong signs. This 
gives people an impression that the coefficients Li, ■ ■ ■ , Liq might mainly be contributed by the anomaly part, 
and the normal part might only contribute small corrections However, we note that the results in Eqs.((^) are 

independent of QCD interactions, i.e., these terms remain unchanged when we switch off the QCD gauge coupling 
constant a s . This is somewhat confusing since these coefficients indicate meson interactions which should be residual 
interactions between quarks and gluons. We shall see in the next section that these terms will actually be completely 
cancelled by the terms independent of the quark self-energy in the normal part contribution, so that they do not really 
appear in the final form of the coefficients. What appear in the coefficients are the remaining terms in the normal 
part contribution which depend on the quark self-energy and hence on the QCD interactions as it should be. Another 
feature of the terms in Eqs.(||) indicating that they should be exactly cancelled and should not appear in the final 
formulae for the coefficients is the divergence of Hi and H2 when taking A — > 00. We know from Ref. that the 
ultraviolet divergences in the 0(p 4 ) chiral Lagrangian coefficients come merely from the meson-loop corrections with 
the 0(p 2 ) interactions. In the 1/N C expansion, the meson-loop corrections belonging to 0(l/N c ) will not take place 
in the large-iV c limit. Therefore, in the large- N c limit, the final expressions for the 0(p A ) coefficients should be finite 
when A — > 00. Now the ultraviolet divergences in Hi and H2 in Eqs.(^) have nothing to do with the meson- loop 
corrections, so that they should be exactly cancelled by other terms and should not appear in the final expressions 
for the 0(p 4 ) coefficients. 
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III. ON THE CONTRIBUTIONS FROM THE NORMAL PART 



In this section we use the same regularization technique as in Sec. II to calculate the normal part contributions to 
the chiral Lagrangian coefficients. We start from the effective action S^ ™ 1 ^ given in Ref. ||, 
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The symbols are defined in Ref. ||. 

In the large- N c limit, the integrations in Eq.(|ll|) can be carried out by the saddle point approximation with the 
saddle point equations 
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in which the 0(1 /N c ) term Tj is neglected. Note that the last term in Eq.(|l7|) actually vanishes due to Eq. (|l5|) . We 
keep it here for showing the relation between the effective action g^"™) anc ] it s stationary conditions Eqs.(|l3|) — (|lE|). 

In the large- iV c limit, = <&n c on the right hand side of Eq.(12). The left hand side of Eq.(|l2]) can be carried 
out from Eq.(|l7|) using Eqs.@-(jl5]). Then the explicit form of Eq.(jlf) in this approximation is 



i[(t? + Jn - nnc)" 1 ]^^, x) = $^(ar, x). 



(18) 



We see that Un c and $q c play the roles of the quark self-energy and the quark propagator, respectively, in the case 
with Jo ^ 0. 

Now we decompose »Sgg 0rm ' ) into a part independent o/IIo c and a part depending on line- The part independent of 
ITn c can be extracted from 5^° rm ' ) by setting ITn c = 0, i.e., 
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Here we have ignored the last term in Eq.(p"7|) which actually vanishes due to Eq. (|l5|) . We show in APPENDIX B 
that the last term in Eq. (|l9|) is actually ^-independent. Therefore, Eq.(|l9|) can be written as 

^norm,n nc =o) = m ( + j Q j + q _ ^dependent terms] . (20) 

Comparing the J^-dependent terms in Eqs.(^J) and (p0|), we see that they are of the same form but with an opposite 
sign. Thus their contributions to the chiral Lagrangian coefficients exactly cancel each other to all orders in the 
momentum expansion. The cancellation in the case of the 0(jp 2 ) coefficient Fq has been described in footnote 1 in 
Sec. II. For the 0(p 4 ) coefficients, we have 
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Thus, in the large- iV c limit, the anomaly part contributed chiral Lagrangian coefficients in Eqs.(||) do not really appear 
in the final results of the chiral Lagrangian coefficients although their values Eqs.(^[) are close to the experimental 
values. The chiral Lagrangian coefficients are actually contributed from the lTn c ^ part of S^™\ 
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This is our first new conclusion in this study. 
The final chiral Lagrangian coefficients are then 
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Since IIf2 C | Ss= q = 0, these 0(p 4 ) chiral Lagrangian coefficients will vanish if we switch off the QCD coupling constant 
g s as it should be. 

IV. CALCULATION OF THE CHIRAL LAGRANGIAN COEFFICIENTS 

We see that to calculate the chiral Lagrangian coefficients from c^ 01 " 1 * 11 "^ ) ^ we should mainly deal with <S^ orm ' 
given in Eq.(^) which has never been carefully calculated in the literature. Ignoring the vanishing last term in 
Eq. (|l7j) , there are still rather complicated terms in it. For example, the third term includes various ranks of gluon 
Green's functions which concern very complicated calculations of QCD dynamics. As the first time of doing this kind 
of calculation, we shall take futher approximations to simplify the evaluation of sjS >TW ^ ■ We know that the pion decay 
constant f n has been studied from QCD in Ref. by taking the approximation of keeping only the leading order in 
dynamical perturbation, i.e., taking into account only the QCD interaction in the Schwinger-Dyson equation leading 
to the nonperturbative solution of chiral symmetry breaking, and neglecting other QCD corrections in positive powers 
of g s (perturbative). This approximation leads to the widely used Pagels-Stokar formula which is reasonable though 
not perfect. In the large- N c limit, f n is just the 0(p 2 ) chiral Lagrangian coefficient Fq given in Eqs.(^). So, as in Ref. 
H , we take the approximation of keeping only the leading order in dynamical perturbation to calculate the chiral 

Lagrangian coefficients from iS^g 01111 ''. In this spirit, we neglect the complicated third term in Eq.(fL7|) which contains 
only positive powers of g s . Furthermore, we see from Eq.(|l4]) that the second term in Eq.(|l7|) is of the same order as 
the third term, so that we neglect the second term in Eq.(|Tj]) as well. With this approximation, S^g ™^ is simplified 

as 

^norm) = _ ^ + ^ _ ^ _ (2g) 

Now the concerned QCD dynamics resides in the II^c term which is related to the quark self-energy [cf. Eq.([l8])]. We 
expect such a simple approximation may also lead to reasonable results of the 0(p 4 ) chiral Lagrangian coefficients 
since reasonable values of the 0(p 4 ) chiral Lagrangian coefficients have been obtained in a model by Holdom ||] 
considering only the quark self-energy contribution with certain phcnomcnological ansatz. We shall see in Sec. V 
that our obtained 0(p A ) chiral Lagrangian coefficients are indeed reasonable. Although this approximation is crude, 
it provides a simple illustration of the main feature of how QCD predicts the chiral Lagrangian coefficients. Further 
improved study beyond this simple approximation taking into account the second and third terms in Eq.(^) is of 
course needed. That will be presented in a later paper. Now we need to calculate D^(x, y) and carry out the explicit 

expression for Sj^ lm ^ in Eq.(|26|). 

We have noticed that is related to the quark self-energy. If we find out the relation between n^(x,j/) and 
the conventional quark self-energy £(— p 2 ), then we can obtain £(— p 2 ) by solving the well-known Schwinger-Dyson 
equation. As in the literature, we shall write down the Schwinger-Dyson equation in the Landau gauge which is stable 
against the gauge parameter. In the same approximation of taking the leading order in dynamical perturbation theory 
and with the improved ladder approximation, the Schwinger-Dyson equation in the Euclidean space-time reads P-[Tl|] 

2 , 3iV c f d 4 q a s [(p - q)] X(q 2 ) 



S(p ) - — y ^3 {p _ q) 2 q 2 + ^2 {q 2 ) - - ( 27 ) 

This equation can be solved numerically, and the details will be presented in Sec. V. Naively, we may expect that 
H t Q c {x,y) = S' Tp ^(d 2 )S 4: (x — y). But this is not correct. Under a local chiral transformation h(x) (hidden symmetry 
transformation ||), transforms as 

n nc {x,y) ^U' nc (x,y) = h)(x)U nc (x,y)h(y), (28) 

while S ap Y}(d 2 )S (x — y) does not transform like this. The correct relation can be found by replacing the ordinary 
derivative d£ by the covariant derivative 

V£ = fl£ (29) 

Since the external source Vq(x) transforms as 

v&x) - <(x) = hHx)v^x)h(x) + ih\xWh(x)], (30) 
the covariant derivative V„ transforms as 
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Thus the correct identification is 



(31) 
(32) 



Then the effective action (feq) can be written as 
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Next we evaluate the effective action ( |33| ) using the Schwinger proper time regularization as before [cf. APPENDIX 
A] to obtain the expressions for the chiral Lagrangian coefficients. This is not trivial since usually this regularization 
scheme is used in the case that £ is a constant, and thus the momentum integration can be explicitly carried out 

2 i 1 

to check the local gauge invariance of the result. Now we leave E(V ) as an unspecified function in Eq.(33), so that 
the momentum integration cannot be carried out explicitly. Organizing terms to guarantee local chiral invariance is 
rather tedious and the details are given in Ref. j^] . Our obtained expressions in the Minkowskian space-time are 



F£B = 4 



*o =2 



^.(norm) 



x„ 



(-2E» - p%£ p )A 2 + (2S 2 +p>E p V p y-jJL 

J X 2 F_ r, 2 

dp 



2A p Xl 



X 1 Y r, 2 Y r. 2 

9A —L-A P I P P y/2 X 2 

^2 p A 4 2 p A 2 2 p P' 



KX 



(norm) 



dp 



dp 



- 2B p X p 
11 3 



A 2 



B^ + C^X^X 



P A 4 



2p ^K P \(6X± 



6X1 
A 2 



-(-4£ 2 + ^)(-2X p 3 



2X 



X„ 



X„ 



P _ ^P_\ _ , y2 

A 2 A 4 A 2 p 



A 



(norm) 



dp 



-4S p 



2p 2 S 2 



6X1 
A 2 



A 2 
3A 2 
A 4 



3A 2 
A 4 



2 -p-p 
Xp 



x 



^)+4£ 2 (-2A 3 



2X 



A 2 



A 4 A 2 



(norm) 



(norm) 



(norm) 



K. 



(norm) 



A 



(norm) 
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jr^(norm) 



= /C, 

= 2 

= 
2 


4 



(norm) 



dp 



dp 



dp 



= 0, 
(3S 2 



X„ 



■ 2p 2 E p £ p )X 2 + [-2E 2 -p 2 (l + 2E P 3£)]- S f 



(S 2 + 2p 2 E p E p )X 2 - p 2 (l + 2£ p £ p ) 



A p 
A 2 " 



(-4£ 3 +p 2 £ p )X 3 



(4£ 3 



P%)^f -(2£ 3 



l n 2 Y + 3y ^£ 



3S p A p 



(norm) 



= o, 



A 



(norm) 



/c 
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A 2 
^A 2 " 

(norm) 



dp 



(oP S p S p 



1 

A p 
* AT 

dp 



2F„A, 



p vi p 



2F ^ 
p A 2 



(C p — ^p)"rf — (^p 



AO*; 



2EpX p 



p A 4 



P W2 *P P v /2 x^2 

2 ^A 2 " ~2 p p 



(34) 
(35) 
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(36) 



in which the short notations (in the Minkowskian space-time) are 

£ p ee £(-p 2 ), (37) 

x 1 



p ~ p 2 - £2 



p 



a p = -- p 2 e p e p (-i - 2e p e;) - -s 2 (-i - 2e p e;) + - P 2 s 2 (-s; 2 - e p e;') 
-i P 4 (-s; 2 -s p s;'), 

5 p - -Hp%£;(-1 - 2S p Ep - iE 2 (-l - 2£ p £ p ) + \p 2 Y? p (-^ - S p Sp 
-^p 4 (-S; 2 - E p ££) - ip 2 (-l - 2£ p £ p ), 



r - i - in v' - ir) 2 r' 2 
p — 3 3 p p 2 p ' 

Dp = l P 2 s; 2 - ^%s;'(-i - 2E pS p - -/^'(-i - 2 s p sp - ^ 4 s; 2 (~s; 2 - s p s;o 
-i P %£ p (-£; 2 -£ p £;')> 

= -^%£ p (-i - 2s p sp 2 - ^ 4 s; 2 (-i - 2£ p £ p ) 2 , 

^ = + ^ 2 (s p £ P ) 2 - + \^' P + ^ 2 s 2 (-s; 2 - s P sp 

™p 4 (-£? - E„££) - ip 2 (-l - 2£ p £ p ) - ip 2 . (38) 

For the coefficient Fq, Eg. (p5|) is just the well-known Pagels-Stokar formula || when taking the regularization cutoff 
parameter A — > oo. 

It is easy to check that these /Q nonn ) (i = 1 ; • • • 15) do contain the ITn c -independent (£ p -independent) piece which 
exactly cancel the anomaly contributions in Eqs.(Q) mentioned in Sec. III. This can be done by taking a constant £ p to 
carry out the momentum integrations, and picking up the E p -independent terms which are just the IIn c -independent 

terms mentioned in Sec. III. Subtracting these IIf2 C -independent terms from the obtained /(;( noim ) m Eqs.(f36|), we get 
the desired ^( norm ' nn <=^ ) m Eq.(p3|), which is needed in the final expressions for the chiral Lagrangian coefficients in 
Eqs.©. 

We can also check that the regularization cutoff A does not appear in /^j* 10 ™ 1 ' 11 "^ ) ] so that the obtained chiral 
Lagrangian coefficients L\, ■ ■■ , Liq are all finite as it should be since there is no divergence in the large- N c limit 
[the divergent meson-loop corrections are of 0(1/N C )]. 



V. NUMERICAL CALCULATIONS 



The last step in the calculation is to solve the Schwinger-Dyson equation ( p7| ) numerically to obtain E(p 2 ). In 
the integrand in the Schwinger-Dyson equation (p7j), there is still the QCD running coupling constant a s (p — q) 
unspecified. The high momentum behavior of a s is well known. The one-loop level formula is 

/ \ p 2 ^oo 127T 1 

ocsip) — ► 7^ — ^7T , -, 2/A2 r< (39) 

(33 - 2N f ) m(p 2 /A£ C£) ) 

The low momentum behavior of a s {p) is not known yet due to the ignorance of nonperturbative QCD. Inevitably, we 
have to take certain QCD motivated model for it as in the literature. We shall take the following Model A from Ref. 
and Model B and Model C from Ref. |l(]] as examples to do the calculation. They are 



9 



A : a s (p)= 7 /oo OArV for hV/A 2 ^) < -2; 



12tt 
(33^2^,) 

{7-|[2 + ln(p 2 /A^ cg )] 2 } (3 3^ } , for - 2 < \n{p 2 /A 2 QCD ) < 0.5; 



1 12tt 



ln(pVA 2 ?CI) )(33-2Ar / )' 



for 0.5 < \n(p 2 /A 2 QCD ). (40) 



B: a.(p)= 47rW* 4 (p) + ,„„ * 2 ?L , W9 - * . 2 r ? (41) 

(33 - 2N f ) ln(2 + p 2 /A 2 QCD ) 

^ , s 47T 3 , „2 /„2 127T 1 , , 

C: Qs P = — P e + ^ ovn 2/A2 V 42 

p (33-2AT / )ln(2 + p 2 /A^ CIJ ) 

They all have the asymptotic behavior (|39|). In Eq.(40), there is only one parameter Aqcd, while in Eqs.(^Tj) and p^), 
in addition to Aqcd, there are extra parameters 77, and respectively. We shall determine the parameters in the 
following way. In the present approach, there are no meson-loop corrections. Thus we should identify Fq = f w = 93 
MeV ||], and Fq is given by the Pagels-Stokar formula. Changing the parameters will cause a change in £(p 2 ), and 
thus a change in Fo. We take Fq = 93 MeV as a requirement to determine the parameters. In the case of Model A, 
the determined Aqcd is Aqcd — 484 MeV ||. In the cases of Model B and Model C, there are extra parameters. 
We take the original values pq — 380 MeV, and Aqcd — 230 MeV as in Ref. p0[ , and determine 77 and p in the above 
way. The determined values are 77 = 290 MeV, and \i = 1160 MeV [J , The running coupling constant a s {p) in the 
three cases are plotted in FIG. 1. We sec that they are different mainly in the low momentum region. 

To solve the Schwinger-Dyson equation ( |27| ) , we further take the usual approximation a s {p~q) ~ 9(p 2 — q 2 )a s (p 2 ) + 
6(q 2 — p 2 )a s (q 2 ) |l2| with which the angular integration can be easily carried out, and the integral equation can be 
converted into the following differential equation ^] 



d P 2 ( c s (p) \ ' &tt P 2 + £ 2 (p 2 ) 
p 



V' i . « = 0. (13) 



E( AV^^/V *5g> =0, (44) 



with boundary conditions: 

Y.(\ 2 ^-- , ,,„ 

7r A 2 Jq y q 2 + T 2 {q 2 ) 

S'(0) + = 0, (45) 

where A is a momentum cutoff regularizing the integral. We shall eventually take A — > 00. 
We know that the asymptotic behavior of S(p 2 ) reflecting chiral symmetry breaking is 

, ln''- 1 (p 2 /A 2 CD ) 
S (p2 ^ ; Q CD \ (46) 



where 7 = (9iV c )/(2(33 - 2iV/)). We have found the numerical solution of Eqs.([43j), Q44J) and (|45j) satisfying this 
asymptotic behavior. The obtained solution with A — > 00 (a large enough number which can be regarded as infinity) 
in the three cases are plotted in FIG. 2. Again they are different mainly in the low momentum region. 



2 The original values of 77 and ji in Ref. |10| are 77 = 920 MeV, fi = 600 MeV which are different from ours. The reason is that 
in Ref. JlCI] the number of quark flavors is taken as Nf = 6 rather than Nf = 3, and the formula for f„ is more complicated 
than the Pagels-Stokar formula. 

3 In the case of Model B, there is a term containing 8 4 (p) which is not a function of p 2 , and the integration can be directly carried 
out. Therefore, in this case, the differential equation and boundary conditions are different from Eqs.(^), (^), and (Eq). They 

s(p 2 ) 



r s( P 2 )[i- 



7] 



dp 2 



dp z p z lnC2 + p 2 /A 2 



3Nc „2 
8tt f 1? 



+£ 2 (p 2 ) 



0, 



E(A J )[1- 



3iV c /2 1 
A 2 +E(A 2 )1 



3N C a a (A 2 ) 
Sti A 2 



f 

JO 



rf„ 2 XH£> 



0, and 



+ 



; : ;-V, 



a s (0) 



E(0) 



= 0., respectively. 
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With the obtained £(p 2 ), we can calculate the 0(p 4 ) chiral Lagrangian coefficients from Eqs.(|5), (p3|), and (|3£|). 
The obtained values of Li, ■ ■ ■ , Liq are listed in TABLE I together with the experimental values |^| for comparison. 
Note that there is no running of L\, ■ ■ ■ , Liq in this simple approach since the meson-loop effects causing the running of 
L\i ■ • ■ , Lio J2| are of the order of 1/N C , and are neglected in this approach. Thus the predicted numbers of Li, ■ ■ ■ , Liq 
can be directly compared with the experimental values. We see from TABLE I that: 

i) these coefficients are not so sensitive to the forms of a s (p)', 

ii) all the obtained L\, ■ ■ ■ , Liq are of the right orders of magnitude and the right signs; 
hi) Li, L2, £4, Lq, and Liq are consistent with the experiments at the la level; 

iv) I/3, L5, Lj and Ls are consistent with the experiments at the 2a level; and 

v) only Lg deviates from the experimental value by (3 — A)a. 

Considering the large theoretical uncertainty in this simple approach, the obtained L%, ■ ■ ■ , L\q are consistent with the 
experiments. We see that the nonperturbative quark self-energy plays an important role in QCD contributions to the 
chiral Lagrangian coefficients. This supports the phenomenological model of Holdom 

In addition to L%, • • • , Liq, we can also calculate the quark condensate ($ip) from the 0(p 2 ) coefficient FqBq in 
Eq. (|34|) . In the simple approach in this paper, the relation between {tpip) and FqBq is || 

<^A) = -NjFqBq. (47) 

We know that, in this simple approach, Fq = — 93 MeV is finite. But FqBq in Eq.(|34"|) is divergent, 

^ 2 B (A 2 ,A^ CD ) cx \n< (48) 

\ 1V QCD/ 



so that it needs to be renormalized. We take a simple renormalization scheme by taking the counter term as 
Fq Bq(A 2 , /i 2 ), in which [i is the renormalization scale 0. Thus the renormalized quantity 

F 2 B 0r cxlxp( T ^-). (49) 

\ QCD / 

Then the renormalized (^) is 

(^V)r = ~N s FIBq t . (50) 

We take the renormalization scale to be fi = 1 GeV to define the quark condensate, the obtained values of ("tp^r for 
the three forms of a s (p) are 



A : 


(i/np) r = 


-(296 MeV) 3 


B : 




-(296 MeV) 3 


C : 


(W)r = 


-(301 MeV) 3 



(51) 

These are to be compared with the experimentally determined value (ipip)expt = —(250 MeV) 3 from the QCD sum 
rule at the scale of the typical hadronic mass |14|. Considering the large theoretical uncertainty in this calculation, 
the predicted quark condensate is also consistent with the experiment. 

The above results show that the present simple approach does reveal the main feature of the QCD predictions for 
the chiral Lagrangian coefficients although the approximations in this approach are rather crude. Of course, further 
improvements of the approximations beyond this simple approach are needed. This kind of study is in progress. 

Finally, we would like to mention that, in our calculation, we have taken the ultraviolet cutoff parameters A, A — > 00, 
i.e., we have taken account of the QCD contributions in the whole momentum range. Note that this has nothing to 
do with the validity range of the chiral Lagrangian determined by the range in which the expansion in the meson 
momentum makes sense, i.e., up to A x w 47r/ w . To see the role of the QCD contributions from the high momentum 



4 This corresponds to the modified minimal subtraction (AIS) scheme |13| 
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region, say above 1 GeV, we have made a check by doing the calculations with the same £(p 2 ) but taking A = 1 
GeV instead of A — > oo. The results are listed in TABLE II. Comparing the nonvanishing results in TABLE II with 
the corresponding A, A — > oo results in TABLE I, we see that this change of A does not cause much difference in 
L5, L7, L$, and Fq, while it causes Li, Li, L3, Lg, and Liq to reduce by at least a factor of 2. Therefore, we see 
that L5, Lj, L%, and Fq are mainly contributed by the QCD dynamics in the low momentum region, while high 
momentum region contributions to Li, L2, L3, Lg, and Liq are not negligible. 



VI. CONCLUSIONS 



In this paper, we have calculated the coefficients in the Gasser-Leutwyler Lagrangian from the underlying theory 
of QCD in a simple approach with the approximations of taking the large- N c limit, the leading order in dynamical 
perturbation theory, and the improved ladder approximation based on the QCD formulae given in Ref. Q to illustrate 
the main feature of how QCD predicts the chiral Lagrangian coefficients. In the calculation, we use the same regular- 
ization technique, the generalized Schwinger proper time regularization, in the calculations of the contributions from 
both the anomaly part and the normal part, so that the relation between the contributions from the two parts can 
be clearly seen. 

We first take the large- N c limit to evaluate the effective action in QCD. Our first conclusion in this study is that, 
in the large-N c limit, to all orders in momentum expansion, the anomaly part contributions to the chiral Lagrangian 

coefficients [cf. Eqs.^)] given in the literature ^ from the effective action 5^ nom ^ [cf. Eq. (j^)/ are exactly cancelled 
by the contributions from the piece of the effective action c^ olm > n ^-0) y c yr Eq.ftn\)] in the normal part contributions, 
so that the chiral Lagrangian coefficients are eventually contributed by the remaining piece of the normal part effective 
action < § , ^ orm ' nf!c ^ ' f c £ Eq.^Sj]. The final QCD expressions for the 0{p A ) coefficients are given in Eqs.(p5|). 

To simplify g^ 01111 * 11 "^ ^ we f ur ther make the approximation of taking the leading order in dynamical perturbation 
theory. Then g^ ™' 11 ^^ ) j g rec j uce( j to the simple form in Eq.(|33|), and all the chiral Lagrangian coefficients are 
approximately expressed in terms of the quark self-energy £(p 2 ) shown in Eqs.(|34|) — (|36|). To solve the Schwinger- 
Dyson equation for £(p 2 ), we further take the improved ladder approximation. Lacking of the knowledge about the 
running coupling constant a s (p) in the nonperturbative region, we take certain models for it from the literature |9||lC|] 
[cf. Eqs.([40|), ( ^l|) and fll^)], and we further take the usual approximation a s (p—q) « 6(p 2 — q 2 )a s (p 2 )+6(q 2 — p 2 )u s {cf) 
to simplify the calculation. The quark self-energy reflecting chiral symmetry breaking is obtained by solving the 
simplified Schwinger-Dyson equation numerically. The obtained results of the 0(p A ) coefficients are listed in TABLE 
I. Compared with the experimental values of L\ ■ ■ ■ , L\q J2j, the agreement of L\, Li, L4, L§, and Liq is of the 
level of la, and that of L3, L5, L7 and L$ is of the level of 2a. Only Lg deviates from the experimental value by 
(3 — 4)er. Considering the large theoretical uncertainty in this simple approach, all the obtained coefficients L\ ■ ■ ■ , L\q 
are consistent with the experiments. We have also calculated the renormalized quark condensate (V"A)r from the 
obtained 0(p 2 ) coefficient [cf. Eq.(|51)] which is also consistent with the experiment. 

Although the approximations in this simple approach are rather crude, the above results show that this simple 
approach does reveal the main feature of QCD predictions for the chiral Lagrangian coefficients. For studying physics 
not requiring high precision, this simple approach may already be useful. Of course further improvements of the 
approximations beyond this simple approach (reflecting more about QCD dynamics) are needed. This kind of study 
is in progress and will be presented in another paper. 

The approach can also be applied to electroweak theories to study how the coefficients in the electroweak chi- 
ral Lagrangian are predicted by various kinds of underlying gauge theories of the electroweak symmetry breaking 
mechanism. This kind of study is also in progress, and will be presented in separate papers. 
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APPENDIX A: FUNCTIONAL DETERMINANT CONTAINING QUARK SELF-ENERGY 



In this appendix, we take the Schwinger proper time regulation to regularize the one-loop functional determinant 
in which the quark self-energy E(V 2 ) reflecting chiral symmetry breaking takes place. 
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For convenience, the evaluation is done in the Euclidean space-time, and will be analytically continued to the 
Minkowskian space-time after the evaluation. The functional determinant is complex. The imaginary part is just 
the Wess-Zumino-Witten term, and its expression in terms of E has already been given in Ref. Jlf| which exactly 
coincides Witten's result ju|. The phenomenology of the Wess-Zumino-Witten term is well-known and is not related 
to the main purpose of this paper. So we shall ignore the imaginary part here, and concentrate on the evaluation of 
the following real part of the functional determinant 



RelnDct[ J D + E(-V 2 )] = -Trln 



[£>t + E(-V ")][D + S(-V )] 



-- lim f°° - Tr e^lE-v'+s^-vVfcEt-v^+st-v^fl-^t-v 2 )] 

2 A— >oo li T 



(Al) 



where 



D = y - sq + ipnj5, 



dp, - i% - £acty75 = -Vj,, = - iv£(x), 



E-V 2 + E 2 (-V ) + Io£(-V 2 ) + E(-V 2 )/o - I £(-V ) = [D* + £(-V )][D + S(-V )], 
In = -i&nl5 - so - ipnl5, In = -ifisfYs - sq + ipnlz, 

[4 E(-V 2 )] = 7"[d,S(-V 2 )] - -f (d^(-V 2 ) - i[vau, £(-V 2 )A , 



The matrix element in Eq.(Al) can be evaluated in the momentum representation 



(x\e 



-r[B-V 2 +S 2 (-Vl+/nE(-V")+S(-V' I )/n-ji S(-V z )] | 



d 4 p 



exp < — t 



—2 



E(x) - V 2 - 2ip • V K + p 2 + E 2 (-V - 2ip • + p 2 ) 



—2 



—2 



InSC-V,. - 2ip • V K + p 2 ) + £(-V x - 2ip • V* + p 2 )/o - ^ E(-V B - 2tp ■ V x +p 2 ) 



(A2) 



(A3) 



Then after lengthy but elementary calculations and expanding in powers of the external sources, we can identify the 
expressions for F$, FqB , ]C^ wlm \ • • ■ , /C^ 011 "' by comparing with the form of Eqs.(|l|), and the obtained results in 
the Minkowskian space-time are just those given in Eqs.(|3q), (pi) and (36) in the text. The details are given in Ref. 



For the evaluation of the effective action S^S 10 ™ 1 ' in Eq.(||) in the Minkowskian space-time, we note that there is 



=2. 



no E(— V ) term in Eq.(Al), but we still have to replace £(— V ) by an infrared cutoff parameter k in Eq.(Al) to 
regularize the infrared divergence. Then the momentum integration can be explicitly carried out, and we obtain the 
results in Eqs.(m) in the text. 



APPENDIX B: SI- INDEPENDENCE OF THE LAST TERM IN EQ.(19) 



Here we show that the last term in Eq. (|9|) , 



n=2 



IInc=0 



with the saddle point equation [cf. Eq.(|l3|)] 



-im +Jn)- l r(y,x) 



(Bl) 



(B2) 



rio c =o 



is ^-independent. The O-rotated quantities in Eqs.(Bl) and (B2) are defined by 



Jn(x) = [n(x)P R + n\x)P L ] [J(x)+ip][Sl(x)P R + tf(x)P L ], 
Qft(x,y) = [tf(x)P R + Sl(x)P L ] $ T (x,y){rt(y)P R + n(y)P L }, 



(B3) 
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and the 2n-point Green's function G^ (xi, x' x , ■ ■ ■ , x n , x' n ) is define by (|] 



"pi—IM 



,fe.-^)K' 1 (a!i)(fU'f^W] • ' ' Kl(^)(^) Q „ /3 „7 M "^:(x„)] 



(B4) 



First we see from Eqs.(|B3l) that Eq.(B2) can be written as 



= (<& T )&*) 



n fic =o 



n nc =o 



[^(^Pfl + n( w )Pi] $r(»,i)[fi t (i)p B + n(i)p L ] 



n c =o 



p(y)P R + n(y)P L ] [(^ +J)- 1 ](i/,a:)[nt(a:)PK + n( a !)P i 



-i per 



in which 



satisfies 



and 



is ^-independent. 



-ip(y)P R + n(y)P L y?'i(i$ + J)- 1 ]?'*' (y,x)[tf(x)P R + n(x)P L f c 

-iboVmo] w (v)[W + J)- 1 / a '(y,x)v^'(x), 
v a (x) = afi(x)p R + n(x)p L 

7oV r s ](x)7o7 Al = l^{x), 

= -i[(i$ +J)- 1 }(y,x) 



(B5) 
(B6) 
(B7) 

(B8) 



n c =o 



With the expression (B5) for <&ff c (x, y)|n s -2 C =o, Eq.(Bl) becomes 

4 , HWcS?)"- 1 



n=2 



n nc =o 



(aii,^,-- • ,a;„,a;^)(7oV^7o)' Tl<Tl (a;i) • ■ 



X (70^70)^^ (z n )V^"Vi) • ' ' V^ p "(x' n )^(x 1 ,x' 1 ) ■ ■ ■ $^(x n ,x' n ) 



n c =o 



n=2 



N C J2 d i x 1 ---d i x' n 



, (-*)"(^v c5s 2 )"- 1 



-G 



(x 1 ,x[ : ■ • ■,x n ,x' n )®Z lP '-(x 1 ,x , 1 ) ■ • ■ $c" P " (x n , x' n ) 



(B9) 



n c =o 



where G^"> , (a?i, a?i, • • • , £„, <) is 



Va,p't~-p'„ 

r Vn,pi—p. 



GS."; ff ™..„ (xi,^, • • ■ ,x n ,x' n ) = (7oV^7o) <TlCTl (a;i) • • ■ (70^70)^" (x n ) 

xG;/;;.;f a^, • • • , Xn , x'jv^ K) • • • v£»> «). 



(BIO) 



Next, we look at this transformed Green's function Gy n p x[, • ■ • , x n , x' n ). From the definition (B4) and the 
property ( |B7| ) we have 



14 



I • • • d 4 <sr 2 G^X". Pii ( Xl , x',, ■ ■ ■ , x n ,x' n yf a \ { Xl )r a \ K) • • • CI Mr a n n K) 

= Gj.V-X^'- • • ^«)[(^7o^7o)S 1 1 ( a; i)(^ 1 )« 1 /3 1 7 Ml (^V')S(xi)] • • ■ 
• ■ • [@l*V^ a )ll{x n )(^) anfinl ^ (Va^) a l(x n )] 

= ■ ■ ■^nM\(xi)(^u 101 r^ i 1 (^)} ■ ■ ■ [Tz{x n ){^) an0nl ^r P :{x n )] 

= J d 4 x[ ■ ■ ■d 4 x' n gl l ~ 2 G a p 1 1 '.'.'.p^{xi 1 x / 1 , ■ ■ ■ , x n , x'Jip^] ' "Cj^lftt 1 !.). 



r^O\ ■■■<T n 

tJ Vh,Pl•••p„ , 



• • • , x n , x n ) — G 'p 1 ...p n {x\i ■ ' ' j x n , Xn ) . (Bll) 
Thus the transformed Green's function Gy^'""^ in Eq.(|BS|) can be replaced by G^ 1 .'.'.'^, and Eq. (^9| ) becomes 



n=2 



s^ = N c y2 d* Xl ---d*x> n 



- G pI---p: (Sl. *!> ■ • ■ . X n, <Wc 1Pl {X 1 ,X' 1 ) ■ ■ ■ ( Xn ,x' n ) 



(B12) 



n c =o 



which is independent of f2. 
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Tables 



TABLE I. The obtained values of the 0(p 4 ) coefficients Li ■ ■ ■ , Linfor Model A [Eq.@], Model B [Eq.©] and Model C 
[Eq.(^)] with A, A — > oo together with the experimental values [Eq.(^)] for comparison. Aqcd is in MeV, and the coefficients 
are in units of 10 -3 . 





Aqcd 


la 


L 2 


L 3 


U 


L B 


Lb 


L 7 


Ls 


L 9 


Lio 


A: 


484 


1.10 


2.20 


-7.82 





1.62 





-0.70 


1.75 


5.07 


-7.06 


B: 


230 


0.921 


1.84 


-6.73 





1.43 





-0.673 


1.64 


3.80 


-6.22 


C: 


230 


0.948 


1.90 


-6.90 





1.29 





-0.632 


1.56 


3.95 


-6.21 


Expt: 




0.9 ±0.3 


1.7 ±0.7 


-4.4 ±2.5 


0±0.5 


2.2 ±0.5 


0±0.3 


-0.4 ±0.15 


1.1 ±0.3 


7.4 ±0.7 


-6.0 ±0.7 



TABLE II. The same as in TABLE I but with A = 1 GeV instead of A -» oo. 





Li 


L 2 


L 3 


L 4 


Ls 


L 6 


L 7 


Ls 


L 9 


Lio 


F 


A: 


0.403 


0.805 


-3.47 





1.47 





-0.792 


1.83 


2.28 


-4.08 


88.7 


B: 


0.281 


0.563 


-2.71 





1.44 





-0.836 


1.83 


1.46 


-3.69 


89.6 


C: 


0.304 


0.608 


-2.86 





1.43 





-0.855 


1.87 


1.56 


-3.64 


89.4 
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Figures 




(p 2 ) 1 ' 2 (GeV) 



FIG.l. a s (p) for Model A [Eq.(|0|)], Model B [Eq.(|4lj)], and Model C [Eq. (|42j)]. The solid, dashed, and dotted lines 
are for Models A, B, and C, respectively. 




(p 2 ) 1 ' 2 (GeV) 



FIG. 2. The obtained E(p 2 ) from the Schwinger-Dyson equation @ with Model A [Eq.©], Model B [Eq.@] and 
Model C [Eq. (|4|)] for the running coupling constant a s . The solid, dashed, and dotted lines are for Models A, B, 
and C, respectively. 
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